Entangled photon pairs from a quantum dot cascade decay: 
the effect of time-reordering 
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Coulomb interactions between confined carriers remove degeneracies in the excitation spectra of 
quantum dots. This provides a which path information in the cascade decay of biexcitons, thus 
spoiling the energy-polarization entanglement of the emitted photon pairs. We theoretically analyze 
a strategy of color coincidence across generation (AG), recently proposed as an alternative to the 
previous, within generation (WG) approach. We simulate the system dynamics and compute the 
correlation functions within the density-matrix formalism. This allows to estimate quantities that are 
accessible by a polarization-tomography experiment, and that enter the expression of the two-photon 
concurrence. We identify the optimum parameters within the AG approach, and the corresponding 
maximum values of the concurrence. 
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I. INTRODUCTION 



In view of their peculiar level structure, semiconductor 
quantum dots (QDs) are considered promising sources on 
demand of entangled photon pairsi^ Although alterna- 
tive strategies have been envisaged, based on the use of 
single-photon sources and postselection^ii the possibility 
of deterministically generating frequency-polarization en- 
tangled photon pairs by a single cascade emission has re- 
cently received strong experimental support There, 
the radiative relaxation of the dot from the lowest biexci- 
ton level generates a two-photon quantum state: \ip p h) = 
{\<j> H ;H,H) + \<j> v ;V,V))/V2, where H and V are the two 
linear polarizations, whereas \4>h) and \4>v) refer to the 
spectral degrees of freedom. In particular, \4>h) (\<frv)) is 
the wavepacket resulting from the sequential emission of 
two photons, with central frequencies lo p and cu r (tu q and 
u> s ) (Fig. [1]). Ideally, u>i = uj p =uo q and u>2 = u r ~lo s (anal- 
ogous equations hold for the relaxation rates): therefore, 
| VyJ = ; H) +|V, V))/y/2®\4>), with |0) = \<t> H ) = |0v). 
In realistic conditions, however, the degree of entangle- 
ment is limited by three main factors£&±£ First, an im- 
perfect system excitation results in a finite probability 
that the system does not undergo a single cascade decay, 
thus emitting more (or less) than the two desired pho- 
tons. Second, the coupling of the confined excitons with 
phonons tends to induce a loss of phase coherence in the 
state of the emitted photons. Third, the presence of an 
excitonic fine-structure splitting tends to make photons 
emitted with orthogonal polarizations distinguishable in 
the spectral domain: Shv = ^ p — ^ q — oj s — lj t ^Q, and 
therefore |(0h|0v)| 7^ 1. This provides a which path in- 
formation, which impedes to rotate the H and V compo- 
nents of \ipph) one into another by linear optics elements 
and to observe interference effects between them. In op- 
timizing the entangled-photon source, part of the recent 
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FIG. 1: (Color online) Level structure of the quantum dot. 
The fine-structure splitting and the biexciton binding energy 
are given respectively by 5hv=E3—E2 and Ab=E2 + E3 — 
E4,. In the WG strategy, the ideal case corresponds to the 
biexciton and exciton emission frequencies being independent 
on the polarization: LOi=iu p —Lo q and UJ2 = u> r = u> s (Shv =0). 
In the AG approach the photon emitted in the B —* Xh 
(B — » Xv) transition matches the color of that emitted by the 
Xv — * G {Xh — » G) decay (As=0). Therefore, wi = lo p — lu s 
(cyan arrows) and L)2=u) q =u} r (red arrows). In order to erase 
the which-path information provided by the emission order of 
the lui and U2 photons, which is opposite in the two paths, 
the photons are delayed by 5 X (x = p, q, r, s). 



effort has been concentrated on the solution of the lat- 
ter problem. Possible strategies include quenching of the 
excitonic fine-structure splitting Shv by means of mag- 
netic field&i or ac-Stark effect^, and spectral filtering of 
the emitted photons,- An ingenious alternative consists 
in engineering the system so as to obtain color coinci- 
dence across generation (AG), rather than within gener- 
ation (WG)i 12 ' 13 i 14 There, the QD spectrum is tuned in 
such a way to have vanishing biexciton binding energy 
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(As = E% + E3 — Ei), so that the four emission fre- 
quencies uj x (x=P,q,r, s) reduce to two: u>i = u> p = uj s 
(cyan arrows in Fig. [1]) and 0J2 = Lo q — lo t (red arrows) . 
In the AG scheme, however, the which-path information 
is provided by the order in which the lo\ and UJ2 photons 
are emitted, which is now opposite for the two polariza- 
tions. In order to make the H and V paths spectrally 
indistinguishable, the photons emitted in the four modes 
X should be spatially separated and time delayed by S x . 
The time reordering implements a unitary transforma- 
tion (U = U H ® \H,H){H,H\ + U V ® |V,V)<V,V| ) of 
the two-photons state, such that: \{<Ph \UlfUv\<t>v)\ > 

\(M<Pv)\M 

Hereafter, we investigate the viability and the limits of 
the AG approach. More specifically, we verify to which 
extent the which path information can be erased by in- 
troducing these frequency and polarization selective de- 
lays. To this aim, we derive analytic expressions for an 
entanglement measure (namely, the concurrence, 7 C) of 
the two-photon state, and derive the expressions of the 
delays 5^ that maximize C, as a function of the emis- 
sion rates T x . This is the same as optimizing the unitary 
quantum erasure U of the which-path information, for 
a given source. In addition, we maximize C(5 X ) with 
respect to T x , thus providing indications for the opti- 
mization of the two-photon source. In semiconductor 
quantum dots, the (relative) values of the exciton and 
biexciton relaxation rates can only be engineered within 
a limited range of values. 16 However, such ranges can be 
potentially extended by coupling the QD to an optical 
microcavity. In the weak-coupling regime, the effect of 
the cavity on the dot dynamics essentially consists in en- 
hancing the photon-emission rates of resonant transitions 
(Purcell effect). Therefore, and in order to allow analytic 
solutions, we don't include the degrees of freedom of the 
cavity explicitly, but rather mimic its effect by enhanc- 
ing r x . We also neglect the effect of pure dephasing and 
imperfect initialization of the QD state {i.e., of realistic 
excitation conditions). In fact, the way in which these 
affect the degree of frequency-polarization entanglement 
is independent on the approach, AG or WG. Detailed 
discussions on these effects, can be found in the litera- 
ture^^ 

The paper is organized as follows. In Sec. [II] we intro- 
duce the density matrix approach we use, and the cor- 
relation functions that enter the calculation of the two- 
photon concurrence. Further details on the method are 
given in the Appendixes. In Sec. IIIII we give the expres- 
sions of the concurrence, and maximize it with respect 
to the relevant parameters. In Sec. [IV] we draw our con- 
clusions. 



II. METHOD 

The time evolution of the dot density matrix, pqd, 
is described by the master equation (fi = 1): pqd = 
i[PQD,H] +J2 x £xPQD, where H = J2k=i E k\ k )( k \ and 



\k) are the QD eigenstates (in the following, we take 
E\ =0). The radiative relaxation processes are accounted 
for by the four superoperators in the Lindblad form: 
£ x Pqd = o- x PQD a x ~ ( a x a xPQD + Pqdv x <J x )/2, with 
X = P,q, t, s. Each of the ladder operators a x corre- 
sponds to one of the optical transitions in the four-level 
system: a v = ^T v |2)(4|, a q = y^\S){4\, a r = y/T\: |1)(2|, 
<7 S = V / T7|l)(3|. The quantum dot is initialized in the 
biexciton state: pqd(0) = |4)(4|; in the absence of mul- 
tiple excitation-relaxation cycles, the cascade-emission 
process from such level results in the generation of two 
photons. The following time evolution of pqd can be 
solved analytically (see Appendix A). Within the AG 
approach, the photons are delayed in time by a quan- 
tity that depends on their energy and polarization (Fig. 
[TJ. The resulting relations between the input mode a x 
and the corresponding output modes a a i read (up to a 
common time delay): 

a>Hi(t) = cr p (t - dp) , a H 2{t) = <r r (t - S r ), 
avi(t) = <J s (t - S s ) , a V2 {t) = <r q (t - 8 q ). (1) 

The quantity of interest here is the degree of en- 
tanglement between the frequency and polarization de- 
grees of freedom of the two-photon state. This can be 
computed from their density matrix p p h which is de- 
rived from the QD dynamics (i.e., from pqd) through 
Eqs. ||T]|. In the following, we refer to the basis 
{\H1,H2),\H1,V2),\VI,H2),\V1,V2)}; here, the first 
(second) mode is identified by the central frequency 
lui = Lu p ((x>2 = u r ), which coincides with lu s {uj q ) in the 
ideal case = 0. Within a single cascade decay, and 
in the absence of non-radiative relaxation channels, the 
matrix elements of p p h correspond to the time integrals 
of second-order correlation functions (see Appendix B): 

(ojl,/32|p p/[ | 7 l,<52) = Jdt' Jdr' G aPlS {t' ,r'). (2) 

Here, G lSa p{t,T) = G% 1& {t, |t|), with ij = 12 for r > 
and ij = 21 for r < 0, whereas 

G%p(t,r>0) = {al l {t)al 3 {t+T)a. a (t+T)a Sl {t))- (3) 

After applying Eqs. (P), the second-order correlation 
functions involving the time-shifted ladder operators a x 
are solved by means of the quantum regression theorem 
(see Appendix A). Experimentally, the matrix elements 
of pph can be accessed within a polarization tomography 
experiment .— 

Given the above master equation and initial condi- 
tions, only few elements of the density matrix do not 
vanish identically (see Appendixes). These are the diag- 
onal elements phh = (H 1, H2\p p h\Hl, H2) and pvv = 
(VI, V2\p p h\Vl, V2), and the off-diagonal one pnv = 
(HI, H2\p p h\ VI, V2). As a consequence, the two-photon 
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density matrix reads: 



p P h = 




P VH 








P HH 



(4) 



The degree of entanglement of the two-photon state 
can be quantified by the concurrence (C), whose value 
ranges from to 1 , going from factorizable to maximally 
entangled states. For the above density matrix, it is eas- 
ily seen that C(p p h) = 2|/?hv|- 



III. RESULTS 

A. Two-photon density matrix 

We start by considering the diagonal matrix elements 
of p p h, namely phh and pvv- The contribution to 
Phh = (HI, H2\p p h\Hl, H2) corresponding to the or- 
dered detection of photons HI and H2 is given by the 
time integral of the correlation function 



U HHHH 



(tV) = (al(t)4(t + T)a r (t + T)a p (t)) 

= r p r r exp[-(r p + r g )t-ryr], (5) 



where t' = t 



dp and r' 



The condition 



that the HI photon is detected before H2 (V > 0) results 
into a lower bound for the delay in the input modes a x : 
t > max(0,(5 r — 5 p ). Analogously, the contribution to 
(VI, V2\p p f l \Vl, V2) corresponding to the ordered detec- 
tion of photons VI and V2, Py V , is given by the time 
integral of the correlation function: 

Gfyyyit' ,T') = (<x\(t) ffj ( t + T ) <7 s (t + t) CT q (t)) 

= r g r s exp[-(r p + r g )t-r s r], (6) 

where t' = t+r+5 s and r' = — T+S g — S s . Here, the photon 
order in the output modes (ay\ and ay%) is inverted with 
respect to that of the corresponding input modes (<r s and 
a q ). This condition results in an upper bound for the 
delay in the emission process: < r<S q — d s . 

After time integration in t and r, the above expressions 
yeld the following coincidence probabilities: 



Phh — r p r r 



dt 



dTe -(r r +T q )t-T r r 



max(0,<5p 



r P + r ? 



^-{exp[-r r xmax(0, S pr )}}, (7a) 

p 1 v 2 y=T q T s dt dre-^+ r ^-^ 
Jo Jo 

{l-exp[-r a xmin(0, *„)]}, (7b) 



Tp + Tg 



with 5 XX > = S x — S x l Analogous expressions apply to 
the case ij = 21. After summing up the contributions 
corresponding to the two cases, the diagonal elements 
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FIG. 2: (Color online) Values of \p H v\ = C M (T)/2 as a func- 
tion of T„ p and T qp (where T xx / = T x /r x / ). Fhe value of 
log(T rp ) has been fixed to —1.2 (a), (b), and 1.2 (c). 



in the two-photon density matrix take the simple form: 

Phh=Phh+Phh =t p/( t p+ t i) and Pvv=pyv+Pvv = 
T q /(Tp + T q ). Even though the two-photon concurrence 

depends on the off-diagonal terms of p p h , upper limits for 

C can already be derived from the diagonal elements. In 

fact, being |pw| 2 < Phh Pvv, it turns out that C < Cq = 

2Tp{ 2 /(l + Tpg), with T XX '=r x /r x '. Such upper limit, 

corresponding to the density matrix p p h of a pure state, 

has an absolute maximum of 1 for T pq = 1. The physical 

interpretation of the above inequality is that, besides the 

erasure of the which-path information, a high degree of 

entanglement in the two-photon state requires a balanced 

branching ratio between the H and V decay paths. 

The relevant off-diagonal matrix element of p p h is given 

by the time integral of the correlation functions (Eq. @): 

G H HVv(t',T') = (al(t)at(t+T)a s (t+S ps )a q (t+T-S gr )) 
= exp[A(t,r)+iB(r)], (8) 



where t' = t + 5 P and t' = r + 5 r 



S p . The real and 



imaginary parts of the exponent in the second line are: 

A{t,T)=A -(T p +T q )t-(T p +T q +T r -r s )T/2, (9) 
B(t) = B -(E a -E 2 -E 3 )t, (10) 

where A = {Tp+T q )5 qr /2~T s (8 qr +5p S )/2 and B = (E 4 - 
E2)5p S — Eib qr . The integration intervals result from 
the requirements that, in Eq. ([8]), all times in the input 
modes be positive (t + r > 8 qr , t > —S ps ), and that the 
biexciton relaxation takes place before the exciton one 
(r < Sp S + 5 qr ) ; otherwise the two-time expectation value 
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on the right-hand side of Eq. ([8]) vanishes identically. 
As a consequence, the phase coherence between the two 
linearly polarized components of the two-photon state 
reads: 

POO rSps+Sqr 

PVH = dt dr ^l+^W. (11) 

J max(0,(5 ps ) J 8 qr —t 

A finite biexciton binding energy would result in an os- 
cillating term e lB ^ T \ and therefore in a suppression of 
\ P hv\ for | A s | > (T p + T q + T r - T„)/2. Within the 
WG strategy, the analogous condition reads: \5hv\ ^ 
(r r + r s )/2. If the resonance condition E4 = E2 + E3 is 
fulfilled, then B(t) = Bo, and the (constant) phase of 
GHHw{ti T ) plays no role. The concurrence that quan- 
tifies the energy-polarization entanglement of the two- 
photon state then corresponds to: 



g _ 2 (II x r x) 1 2 [ r -[T s S ps + (T v +T q +T s )S qr }/2, 

CDE l° e + 



where C=-{T p +T q +T r -T ,) /2, E=-(r p +r q -T r +r s )/2, 
and D = T p + T q . We note that C only depends on 8 ps 
and <5 gr , and not on the four delays 5 X independently. 
This is consistent with the intuition that the degree of 
entanglement depends on the extent to which the delays 
make the HI (H2) mode indistinguishable from VI (V2) 
in the time domain. 



B. Parameter optimization 



Given the analytic expression of the concurrence, 
C(5 ps , 5 qr ; r), we first maximize it with respect to the 
delays, as a function of the relaxation rates. The opti- 
mum values of the delays are denoted with S^(T) and 
8™{T), being T = (T p ,T q ,T r ,T s ). The corresponding 
concurrence is 



c M (r) = c[6%(r),6%(ry,T]. 



(13) 



In a second step, we maximize C M (T) with respect to the 
relaxation rate, thus identifying the absolute maximum 
ofC, namely C opt = C M {T opt ). 

The values of the relative delays that maximize C sat- 



isfy the conditions: dC/d5 ps 
expressions read: 



and.dC/d8qr = 0. Their 



C(r) 



C(r) 



j 1 * 



2r, 



In 



2T r 



In 



2r s 



2T r 



(14a) 
(14b) 



In the general case, one can substitute the above equa- 
tions in Eq. (|12p. thus obtaining C M (T). In order to 
maximize the concurrence, we look for the values of T v 




iog 10 (Y M .) 



FIG. 3: (Color online) The solid black lines give the values 
of the concurrence for optimized delays, C M , as a function of 
the ratio T x x i = T x /Y x i and for different constraints, (a) 

r s = r p and r r = r„ with ( x ,x) = (q,p); (b) r s = r q 

and r r = r p , with (x,x) = (P,q)~, ( c ) r <7 = Tp and r s = 
r r , with (XiX') — (Pi r )- The dashed lines are the values of 
the upper limits for the concurrence Co, corresponding to a 
complete cancelation of the which-path information. Photon 
emission rates that are set equal are denoted by arrows of 
equal thickness in the level schemes on the right-hand side of 
each plot. 



that satisfy the equations: dC M /dT x = 0. No value of T 
simultaneously fulfils these conditions. Therefore, we nu- 
merically compute C M (r) for a wide range of relaxation- 
rates values: T rn < T gp ,T rp ,T sp < Tj;. The global 
maximum of the concurrence corresponds to T p = T q and 
r r =r s , with the latter rates much larger than the former 
ones: 



opt 



lim C 

T or .->o 



M 



(T) ~ 0.736. 



(15) 



In order to gain some further understanding, we plot the 
dependence of C on T q and r s , for three different values of 
r r , all in units of T p (see Fig. [5]). On average, C increases 
with increasing values of r r /r p [from panel (a) to (c)]. 
For T r < T p [panel (a)] , the maximum is localized close to 
the point T sp = T qp — 1. The optimum delays in this case 
[log(T rp ) = —1.2] are no longer identical. In particular, 
they are given by: 8™ ~ 0.80/r p and 8™ ~ 2.2/T p . If the 
Xh exciton relaxation rate is larger than the biexciton 
one [see panel (c), where log(T rp ) = — 1.2], the maximum 
of C is localized close to the point T sp — 1.2, T qp = 1. 
The corresponding delays are reduced to: <5*f = Sgf ~ 
0.061/r p . The above examples show that the best values 
of the delays strongly depend on the relaxation rates. As 
to the dependence of C on T, the above plots suggest 
that the two biexciton (exciton) relaxation rates should 
coincide (are correlated). 

As already mentioned, the exciton and biexciton re- 
laxation rates can, to some extent, be engineered in the 
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growth process^* A further tuning can be achieved by 
coupling the dot with an optical microcavity, through the 
Purcell effect. 7 In the following, we further consider the 
dependence of the concurrence optimized with respect to 
the delays, C M (T), after introducing specific relations be- 
tween the parameters T x . These will either be realizable 
by suitable dot-cavity couplings, starting from a situa- 
tion where r x ~ T [cases (i) and (ii)], or correspond to 
a region of specific relevance for the maximization of C 
{Hi). 

(i) In the first case, the two decay processes with equal 
frequencies share the same value of the emission rates: 
r p = r s and T q — T r . This condition might be fulfilled 
by coupling the QD with a suitable microcavity. In par- 
ticular, the MC should possess a mode doubly degener- 
ate with respect to polarization, and in resonance with 
two QD transitions (e.g., p and s), while sufficiently off- 
resonance with the remaining two. The Purcell effect 
would then result in an effective, frequency-selective en- 
hancement of the emission rates T x . The optimum values 
of the delays reduce to: 5™ = (1/T q ) ln(r 9 /r p + 1) and 

$qr = (!/ r p) ln ( r p/ r <z + !)■ After substituting these ex- 
pressions in Eq. (fl2|) . we obtain [black curve, Fig. [3] 
(a)]: 

c (r|x ps , T, r =l) = - - T ^ 2+Tqp/2+1/{2Tqp) ■ (16) 

The extrema of such function are identified by the zeros of 
its derivative with respect to T 9P . The constrained max- 
imum of the concurrence is: C M (r|Y ps , Y qr = l), which 
is C opt (r|T ps ,T ps = l) = l/2, for T qp = 1/2. This value 
corresponds to half of upper limit for the concurrence, 
Co, in the case phh=Pvv=1/2 (gray curve). 

(ii) In the second case, the relaxation rates depend 
only on polarization: r p = IY and r g = r s . Such situa- 
tion can be induced by a cavity with a linearly polarized 
mode, sufficiently broadened in frequency so as to couple 
to both the QD transitions of a given linear polariza- 
tion (H or V) , while remaining uncoupled with the other 
one. Given these constraints, the delays become: 8^ = 

(i/r p )in(r p /r g + 1) and = (i/r g )in(iyr p + 1). 

This results in an expression for the concurrence that co- 
incides with that of above case (i) [Fig. [3] (b)]. In fact, 
from Eqs. |(T2J), (fHaj) . (fT4b)) one can see that Cj^-^ is in- 
variant under the simultaneous exchange T p <— > T g , and 

r r ^r s : c A/ (r g ,r p ,r s ,r r ) = c M (r p ,r ? ,r r ,r s ) . We 

incidentally note that this is not true in general for ar- 
bitrary values of the delays, i.e. if (5 ps , S qr ) ^ (<5 p f , 5 q l). 
For T pq 7^ 1, the upper limit Cq < 1 (gray curve): thus, 
the two-photon concurrence cannot attain its maximum 
value even for a complete cancelation of the which-path 
information, simply because of the asymmetric branch- 
ing ratio between the H and V paths. The difference 
between C and Co, instead, can be ascribed to the distin- 
guishability between the wavepackets relevant to the two 
polarizations. Therefore neither the energy [panel (a)], 
nor the polarization-selective tuning [panel (b)] of the 



QD photon-emission rates allow to achieve high values 
of the concurrence, and specifically to cancel the which- 
path information required in the AG scheme. 

(Hi) In the third case, the biexciton and the exci- 
ton relaxation rates are independent on the polarization: 
T p = T q and r r = r s . The optimum delays then read: 
K I = S Ps = (V r p) ln(l+r p /r r ), resulting in (black curve) 

C M (r|T p9 ,T rs =l) = 2(T pr + l)-( 1+1 / T -). (17) 

The above expression is a decreasing function of T pr ; 
it tends to Cm = 2/e for T pr — * 0, i.e., in the limit of 
biexciton relaxation much slower than the exciton one. 
This limiting value coincides with the global maximum 
that we find for the unconstrained case. Therefore, as 
already reported above and in Fig. [51 the most favor- 
able region in the parameter space T corresponds to the 
biexciton and exciton relaxation rates being independent 
on polarization, with the former ones much smaller than 
the latter ones. Unfortunately, the present case seems to 
be the least feasible. In fact, within the AG approach 
(where uj p = cu s uj q = u> r ), the transitions B — > Xjj/v 
cannot be resolved from the Xjj/y — > G ones, neither 
spectrally nor through polarization. This impedes to op- 
timize the relaxation rates through the Purcell effect in- 
duced by dot-cavity coupling, and forces to rely on the 
engineering of the QD oscillator strengths alone. We fi- 
nally note that the case (r p T r ) = (r g ,r s ) is the one 
considered throughout Ref. 14. There, analogous con- 
clusions are drawn with respect to the dependence of the 
two-photon entanglement on the ratio T p /r r . However, 
we find that the optimized delays differ from those sug- 
gested by Avron and coworkers, apart from the limiting 
case T r > T p , where 5^ = 6^ -> 1/T r . 

IV. CONCLUSIONS 

We have theoretically investigated the generation of 
energy-polarization entanglement of two-photons emit- 
ted by the cascade decay QD within the AG approach. 
As in the case of the WG scheme, the two-photon en- 
tanglement is limited by dephasing and imperfect dot 
initialization. Unlike that case, the concurrence C is also 
limited by the opposite emission order of the u)\ and uj 2 
photons along the H and V paths (see Fig. [T]). A unitary 
erasure of the which-path information can be performed 
by time reordering; 1 ™ Here, we have analytically com- 
puted C as a function of the parameter that determine 
the time- reordering (i.e., the delays S x ) and characterize 
the two-photon source (i.e., the photon-emission rates 
T x ). We have maximized C with respect to 5 X , as a func- 
tion of T x , thus optimizing the erasure process for an 
arbitrary source. The optimized concurrence C(T) has 
then been maximized with respect to the emission rates, 
thus providing indications on the desirable source engi- 
neering. We find that both the energy and polarization 
selective enhancement of the emission rates, that could 
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be induced by suitably coupling the QD with an optical 
microcavity (Purcell effect), are of limited usefulness. In 
fact, the maximum value C = 0.5 corresponds to identical 
rates (r x = r). On the other hand, the absolute maxi- 
mum of the concurrence, C — 2/e ~ 0.736, corresponds 
to biexciton and exciton relaxation rates independent on 
polarization, with the former ones much smaller than the 
latter ones (T p = T q <C T r = T s ). However, within 
the AG approach (where uj p =w a ^=LO q =u r ), the transi- 



tions B^X 



H/V 



cannot be resolved from the X 



H/V 



•G 



ones, neither spectrally nor through polarization. This 
impedes to access the above regime through the Purcell 
effect. 
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APPENDIX A: DYNAMICS AND 
QUANTUM-REGRESSION THEOREM 

Given the initial conditions pqd(0) = |4)(4|, the 
time evolution of the QD density matrix induced by 
the Liouvillian Pqd (t) = i[pQD,H] + J2 x ^xPQD = 
C QD (t)pQ D (0) is the following: 



(4\PQd\4) 
(3\pqd\3) 

P\PQd\2) 



-(r p +r,)t 



r — r 



,-r.t 



,-rvt 



-(r ? +r a )t 



3 -(r P +r,)t 



(1\ P qd\1) = l-J2( k \PQD\k). 



(18) 



fc=2 



The Liouvillian Cqd doesn't couple the diagonal terms of 
Pqd with the off-diagonal ones. Therefore, for the above 
initial conditions, these are identically zero. 

In order to compute the two-time expectation values 
G l i aaa (of = H, V) , we apply the quantum regression 
theorem (see, e.g., Ref. ? ). This states that if, for some 
operator O, the time dependence of the expectation value 
is given by 



(0(t + r)) = ^2a j (r)(O j (t)), 



(19) 



then 



(A(t)0(t + r)B(t))=J2 a > M (*)*(*)> • ( 2 °) 



In the case of G^ aaa , after performing the substitutions 
reported in Eqs. (TT]), the above operators are: O = cr x a x , 

A = (jyi , B = <r x ', with (XtX') — ( r >p) f° r a — H and 
(X,X') = (s,q) for a = V. Besides, O j = \k(j))(l(j)\ 



(with k,l = 1,2,3,4 denoting the QD state); their ex- 
pectation values (Oj) = pi.k correspond to the elements 
of the quantum-dot density matrix. 

In Eq. (fl"9"j) . the expectation values (Oj(t)) give the 
initial conditions of dot state. Therefore, the two- 
time expectation value in Eq. (|2"0"1) corresponds to the 
single-time expectation value of O, for initial conditions 
P'mjt) = m)Oi{t)B{t)) = (m\A(t)\k(j)){l(j)\B(t)\n). 
This provides an intuitive explanation of why most of 
the matrix elements of p p h vanish identically. The ele- 
ment (H 1, V2\p p h\Hl, V2), for example, corresponds to 
the expectation value of cr^Oq — |4)(4|, for a system ini- 
tialized in the state p'q D = a p pQD(t)aj, cx |2)(2|. Such 
expectation value is zero at any r for Cqd(t)\2)(2\: if the 
dot state is initialized to the exciton level Xh at time t, 
it will never be found in the biexciton state at any time 

t + T. 

The procedure is the analogous for the calculation of 
the correlation functions GHHw{t^ T )- After perform- 
ing the substitutions reported in Eqs. ((T|), however, 
this becomes a four-time expectation value: G l ^ HVV — 
(al(tx)(rl(t2)a s {t3)a q (U)) (where, e.g., h < t 4 < t 2 < 
£3). This is computed by applying three times the quan- 
tum regression theorem, with O = u s 



=4, A<V=I, A^ 



B (i) = j, b(2) =< 



(/ is the identity operator). 



9' 



B&=I 



(21) 



APPENDIX B: 



TWO-PHOTON DENSITY 
MATRIX 



In the present conditions, where the QD is initialized 
into the biexciton state |4) and undergoes a single cas- 
cade decay, the probability that of detecting a photon in 
the al mode and a photon in the a2 mode is given by 
the time integrals of the second-order correlation func- 
tions G^ aaa (t, t) . The identification of such integrals 
with the diagonal elements of the two-photon density ma- 
trix, (al, a2|p p /j|al, a2), results from the fact these have 
the same physical interpretation. 

For the off-diagonal terms the validity of Eq. ^ is 
less intuitive. Such corresponds results from the two fol- 
lowing points, (i) The correlation functions G^ g(t,r) 
and the two-photon matrix elements (7I, (52|p p ^|al, (32) 
transform according to the same equations under the 
change of polarization basis in the 1 and 2 modes. 
(m) The off-diagonal correlation-functions, such as 
Gjj HVV (t,T), can be expressed as linear combinations 
of diagonal ones, namely Gj^g a (t, r), where a and (3 
vary over 4 independent photon polarizations (including 
H and V). These are, in fact, the relations that are 
exploited in polarization quantum tomography. 7 There- 



7 



fore, where Ghhvv = J G'hhvv (*> T ) ^ ^ T an< ^ A Hy being 

the transformation changing the polarization basis. 

Ghhvv = £ A «// G^ a (t,r) dtdr 



a/3 



= ^A^(al,/32|p p ,|al,/32) 

a/3 

= (Vl,V2\p ph \Hl,H2). (22) 
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